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ABSTRACT 

We present three dimensional spectroscopy of eleven E+A galaxies at z=0.06- 
0.12. These galaxies were selected for their strong H<5 absorption but weak (or non- 
existent) [On]A3727 and Ha emission. This selection suggests that a recent burst of 
star-formation was triggered but subsequently abruptly ended. We probe the spatial 
and spectral properties of both the young (<1 Gyr) and old (>few Gyr) stellar pop- 
ulations. Using the H<5 equivalent widths we estimate that the burst masses must 
have been at least 10% by mass (Mb urs t ^10 10 M Q ), which is also consistent with the 
star-formation history inferred from the broad-band SEDs. On average the A-stars 
cover ~ 33% of the galaxy image, extending over 2-15 kpc 2 , indicating that the char- 
acteristic E+A signature is a property of the galaxy as a whole and not due to a 
heterogeneous mixture of populations. In approximately half of the sample, we find 
that the A-stars, nebular emission, and continuum emission are not co- located, sug- 
gesting that the newest stars are forming in a different place than those that formed 
<1 Gyr ago, and that recent star-formation has occurred in regions distinct from the 
oldest stellar populations. At least ten of the galaxies (91%) have dynamics that class 
them as "fast rotators" with magnitudes, v/a, \r and B/T comparable to local, rep- 
resentative ellipticals and SO's. We also find a correlation between the spatial extent of 
the A-stars and the dynamical state of the galaxy such that the fastest rotators tend to 
have the most compact A-star populations, providing new constraints on models that 
aim to explain the transformation of later type galaxies into early types. Finally, we 
show that there are no obvious differences between the line extents and kinematics of 
E+A galaxies detected in the radio (AGN) compared to non-radio sources, suggesting 
that AGN feedback does not play a dramatic role in defining their properties, and/or 
that its effects are short. 
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1 INTRODUCTION 

Galaxies with strong Balmer absorption lines in their spec- 
tra, but weak nebular emission (such as [On]A3727A), repre- 
sent a short-liv ed but potentially imp ortant phase in galaxy 
evolution (e.g. iTran et al.ll2003l . l2004l ). These "E+A" galax- 
ies have strong absorption lines (such as H<5), representing 
a stellar population dominated by A-stars, which are ei- 
ther absent or overwhelmed by the much brighter OB stars 



in most galaxies. These signatures suggest that the star- 
formation within the galaxy abruptly ended £l Gyr ago, 
possibly following a starburst phase. It is likely that a va- 
riety of ph ysical mechanisms lead to such a stellar popu- 
lation (e.g. IPressler fc Gurmlll982l; ICouch fc Sharpleslll987l; 



Zabludoff et all Il996l : iPoggianti et al l Il999l ; iBalogh et all 
19991) . but most invoke a major transformation from one 
galaxy type to another, possibly representing an evolution- 
ary link between gas-rich, star-forming galaxies and quies- 
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cent spheroids. Although such galaxies are very rare in the 
local Universe, their short lifetime means they could poten- 
tially represent an i mportant phase in the evolution of mos t 
normal galaxies (e.g. lYan et alJliooll IZabludoff et al.lll996h . 

In order to trace the route by which star-forming galax- 
ies evolve into quiescent systems, there are several issues 
which must be addressed. First, from where in the galaxy 
does the unusual spectrum originate? Is the strong H5 a re- 
sult of galaxy wide A-star populations, or is it due to a het- 
erogeneous mixture of populations, such as a nuclear star- 
burst with a very high equivalent width (>10A) on top of 
an otherwise normal galaxy? Second, what star formation 
history leads to the spectral characteristics of these galax- 
ies? Is the E+A phase related to a significant enhancement 
of the star formation, or simply the rapid decline of star for- 
mation in normal, gas-rich galaxies. Finally, what triggers 
the recent change in star formation history? 

Several physical mechanisms have been suggested for 
the triggering (and subsequent rapid e nd) of starbursts in 
field E+As, including major mergers (|Mihos fc Hernquistl 
1 19961) and unequal mass mergers and tidal interactions 
l|Bekki et al.l 120051 : iPracv et ail 120051 ). Recently, large-scale 
energetic outflows from AGN have been invoked as a route 
to explain many of the properties of local massive galax- 
ies. If post-starbursts are late-time mergers they may have 
passed through their quasar phase and, as a result, the AGN 
may provi de sufficient energy to inhibit and terminate star- 
formation l|Bower et al.ll2006l ). Indeed, post-starburst galax- 
ies have been sugg ested as ideal laborator ies for testing AGN 
feedback models (|Tremonti et al.l 120071 ) , and some obser- 
vatio n s indicate t hat su ch a link exists l|Georgakakis et ail 
120081 : lYan et all [2006a). Finally, in dense environments 
(such as clusters) galaxy harassment and/or interaction 
with the hot intra-cluster gas is also likel y a driving factor 
i|Dressler fc Gunnlfl98^ ; IPracv et al.ll2010l ). 

While high resolution imaging has demonstrated that 
most bright, nearby E+A galaxies are s pheroidal, often wit h 
signs of interaction (|Yang et al.l 12004 iBalogh et al.ll2005l ). 
their dynamics provide b etter constraints on their likely ori- 
gins. [Norton~e£al. (2001) obtained longslit spectra of twenty 
E+A galaxies and used the ratio of the rotational veloc- 
ity to line-width to conclude that most of the galaxies are 
in the process of transforming from rotationally-supported, 
gas-rich galaxies to pressure-supported, gas-poor galaxies. 
However, longslit spectroscopy mixes spatial and spectral 
resolution, and a better understanding of the dynamics can 
be obtained from integral field spectroscopy, where the spa- 
tial and spectral information can be cleanly decoupled. This 
technique also allows the location, spatial extent and sepa- 
ration of any gas, old (K) stars and (the short-lived) massive 
A-stars to be compared. 

An examp le of s uch measurements were made by 
ISwinbank et~all l|2005l ) using the Gemini/GMOS-South 
IFU to map the properties of a H<5 strong galaxy 
(SDSS J1013+0116) selected from the Sloan Digital Sky Sur- 
vey (SDSS). In this galaxy, the velocity field and spatial 
extent of the distribution of H<5 equivalent widths suggest 
that the post-starburst phase arose from a merger between 
a gas-rich spiral and a gas-poor, passive secondary. Cru- 
cially, within SDSS J1013+0116, the A-stars are distributed 
throughout the galaxy, arguing that the starburst occurs in 
the disk of the galaxy during the interaction. The studies 



Table 1: Target Information 
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J0835+4239 


08 35 21.23 


+42 39 36.0 


0.0920 


10.8 


J0906+5221* 


09 06 19.93 
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0.0989 


10.8 


J0938+0001 


09 38 42.91 


+00 0149.0 
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18.0 


J0948+0230 
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16.2 


J1013+0116 
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0.1055 


10.8 


J1242+0237* 


12 42 52.96 


+02 37 00.9 


0.0846 


10.4 


J1348+0204 


13 48 02.21 


+02 04 05.6 


0.0678 


10.8 


J1627+4800 


16 27 55.98 


+48 00 51.1 


0.1254 


13.0 


J1642+4153 


16 42 55.22 


+41 53 35.5 


0.0722 


11.7 


.11715+5822 


1715 46.24 


+58 22 55.4 


0.1271 


15.6 


J2307+1525 


23 07 43.40 


+ 15 25 59.3 


0.0699 


14.4 



Table 1. Co-ordinates, redshifts and exposure times of the galax- 
ies in our sample. The two E+A galaxies labeled with a * 
(J0906+5221 and J1242+0237) are both in cluster environments. 
Note that J2013+0116 was selected to have weak emission lines, 
so is an e(a), rather than E+A, galaxy. 



of iGoto et ail (|2008l ) and IPracv et all l|2009l ) increased the 
number of E+A galaxies with resolved spectroscopy to ~10, 
and showed that wide-spread A-star populations are a com- 
mon f eature, with A-star populations extended on >>kpc 
scales, [pracv et ail (2009) also showed that the majority of 
their E+A galaxies have significant angular momentum per 
unit mass (so-called 'fast-rotators'), casting doubt on equal- 
mass, major mergers on the mechanism which results in 
E+A's since these should leave dispersion-dominated kine- 
matics. 

In order to gain a better understanding of the interac- 
tion between star-formation and gas dynamics within the 
ISM of E+A galaxies, we have obtained three dimensional 
spectroscopy of eleven galaxies, selected from the SDSS in 
the redshift range 0.07 < z < 0.12. The three dimensional 
data are used to map the distribution and kinematics of 
A-stars. We investigate the spatial distribution of residual 
star formation through the [On]A3727 emission and the dis- 
tribution of the young A-stars through the much stronger 
H<5 absorption line. In §2 we present the data reduction and 
analysis. The results are presented in §3. Finally we summa- 
rize our results and present the implications in §4. We adopt 
a cosmology with f2 m =0.27, A = 1 — f2 m and a Hubble con- 
stant of 72 kms~ 1 Mpc~ 1 . In this cosmology, at z = 0.09 (the 
median redshift of our sample), lkpc subtends 0.6" on the 
sky. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Sample Selection 

Our s ample is selected from the catalog of iBaloeh et al.l 
(2005), which is base d on if-band i magin g of 222 nearby H<5- 
strong galaxies from IGoto et alj (120031 ) . This catalog con- 
tains all galaxies for which the rest frame equivalent width of 
the H<5 absorption line is W (US) >4A and W o ([On])<10A 
(with >2cr confidence), as measured from the spectrum us- 
ing a Gaussian profile fitting technique. Of these, we se- 
lect galaxies which also have weak or absent Ha emission, 
W / D (Ha <5A). This helps to exclude galaxies with weak or 
dusty star forma tion, as well as co ntribution from an active 
galactic nucleus (|Yan et al.ll2006bT ). 

Ten E+A galaxies from this sample were selected for 
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Figure 1. The ftf-band velocity— line-width relation of the E+A 
gala xies in our sample co mpared to SDSS early type galaxies 
from iBernardi et alj j2005h , with i\-band magnitudes taken from 
2 MASS. This figure shows that the E+A galaxies in our study 
span ~80-180kms -1 and ^2 mags in K. Our study focuses on 
the brighter E+A galaxies at a fixed velocity dispersion, due to 
the requirement of high signal-to-noise ratios to carry out this 
analysis. 

follow-up IFU observations. These galaxies were selected 
to span different morphological types and environments. 
The former was quantified by fractional bulge luminosities 
(B /T) determin e d fro m the bulge-to-disk decompositions 
of balogh etaH (|2005T ), while the latter was obtained by 
cross -correlating with the C4 cluster catalogue (Mill er et al.l 
2005). In practice, given our high mass cut, most of the 
ca ndidates were bulge -dominated , as shown in Figure 5 
of iBalogh et al. (2005); we only found one clearly disk- 
dominated galaxy (J1048+0230, with B/T=0. 29). Further- 
more, the vast majority of the E+A galaxies in lBalogh et al.l 
(2005) are not associated with clusters; we found only 
two galaxies satisfying our other selection criteria that are 
near clusters. I n add ition, we include J1013+0116 from 
ISwinbank et all [|2005l ). to bring the final sample to eleven 
galaxies, as shown in Table 1. J1013+0116 was chosen for our 
pilot programme, for which we required reasonably strong 
line emission to ensure reliable dynamical measurements 
could be made. Thus, it has the strongest [On] emission 
m our sample (W o [Oii1=5.2 ±0.3A ), and is more properly 
classified as an e(a) galaxy (|Poggianti et al.lll999l ). This is 
therefore a fairly heterogeneous sample, by design. 

In Fig. [T] we show where these galaxies lie in the 
luminosity-line-width p lane of the Faber- Jackson relation 
l|Faber fc Jackson! 1 197ft ). We measure the velocity disper- 
sions from the SDSS spectra, using the cross correlation 
technique described in § 12.21 We use the line width of the Mgl 
triplet as a proxy for the stellar kinematics of an old stellar 
population. The E+A galaxies in our sample span a range 
of velocity dispersions from a ~80-180 kms -1 , and nearly 
two magnitudes in the AT-band. For compari son we also plot 
the e arly-type galaxies selected from SDSS (|Bernardi et al.l 
120051 '). cross matched to 2 MASS. To ensure a fair compari- 
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Figure 2. The observed (u — g) versus (r — i) distribution (cor- 
rected for Galactic reddening), of galaxies from the SDSS (model 
magnitudes). Small points represent a random subset of galaxies 
from DR5 with 16 < r < 17 in this redshift range; the bimodal 
nature of the population is evident. Filled, coloured circles repre- 
sent our GMOS IFU targets; the red points are t he E+A galaxies 
and th e blue point is the single e(a) g alaxy fromlSwinbank et al.l 
i2005l ). The lines are the predict ions oflBruzual fc Charlotl J2003| V 
models, similar to those used in lBalogh et al. (2005). All assume 
a Salpeter IMF and represent observed colours at z = 0.09. The 
filled, black circle represents an old population, resulting from a 
single burst of star formation 13.7 Gyr ago, while the star in- 
dicates the colour of a galaxy that has undergone constant star 
formation for the same amount of time, and includes internal ex- 
tinction of T v = 1 mag. The solid line represents the track of an 
exponentially-declining star formation rate, with r = 4Gyr, and 
t v = lmag extinction. The dotted line is a "truncated" model, 
which shows the colour evolution over a 2 Gyr period following 
truncation of star formation in the model represented by the open 
star. Finally, the long-dashed line represents a model in which a 
10 per cent (by mass), instantaneous burst is superposed upon the 
old stellar population (filled circle) ; the evolution from blue to red 
is followed for 2 Gyr following the burst. This latter model (which 
includes t v = 1 mag extinction) provides a reasonable match t o 
the colours of our E+A galaxies, as shown in Balogh et al.1 l l2005t) . 

son with the E+A galaxies, we limit the redshift range of the 
comparison sample to z < 0.1, and apply no k-corrections 
to either sample. As Fig. [1] shows, the E+A galaxies in our 
sample are systematically brighter in the 7^-band for a fixed 
velocity dispersion (AMk=0.6±0.2). A compa rable offset in 
_ R-ban d (AniR ~0.6 mags) was also noted bv lNorton et al.l 
|200ll) who used longslit spectros copy of twenty E+A ga lax- 
ies from the sample described bv lZabludoff et all (|l996h . 

The (u — g) versus (r — z) colours of the galaxies in our 
sample are shown in Fig. [2] We compare these to the galax- 
ies from SDSS DR5 with redshift range 0.05 < z < 0.1 and 
magnitude range 16 < r < 17. We also compute three evolu- 
tionary colour tracks: i) an exponentially-declining star for- 
mation rate, with r = 4 Gyr, and t v — lmag extinction; ii) 
a "truncated" model which shows the colour evolution over 
a 2 Gyr period following truncation of star formation in the 
model; iii) a model in which a 10 per cent (by mass), instan- 
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taneous burst is superposed upon the old stellar population. 
As Fig.[2]shows, this latter model provides the best match to 
the co lours of our E+A galaxies, as shown in Balogh et al.l 
(l2005h . 

Finally, we note that of ten galaxies in our sample 
that are covered as part of the VLA Faint Images of the 
Radio Sky at Twenty-Centimeters (FIRST) survey, two 
are detected at >5sigma (J0948+0230 and J1013+0116) 
and two more are detected at 3-5cr (0835+4239 and 
J1642+4153). These detections have 1.4 GHz fluxes in the 
range 0.5-2.2 mjy, corresponding to luminosities of Li.4=8- 
15xl0 21 WHz" 1 . We return to this in 54. 



2.2 GMOS Spectroscopic Imaging 

Spectro-imaging observations of ten new E+A galaxies in 
our sample were taken with the GM OS-North and GMOS- 
Sout h IFU between 2005 and 2007 (|Allington-Smith et all 
2002). Seven targets were observed in 'stare' mode using 
one-slit mode which results in a field of view of 5" x7", while 
four targets were observed with nod-and-shufflc in two-slit 
mode (resulting in a field of 5" x5"). For observations taken 
with nod-and-shuffle we ch opped away from the target by 
30" every 30 seconds (see ISwinbank et al.l 120051 for a de- 
tailed discussion of the observing procedure). The exposure 
times for each target were typically 11-18 ks, with longer ex- 
posures taken on fainter targets (Table 1). All observations 
were carried out in dark time, typically with 0.6" seeing in 
1^-band and in photometric conditions. The corresponding 
physical resolution is ~1.0kpc at the median redshift of our 
sample. We used the B-band filter in conjunction with the 
B600 or B1200 grating, which results in a spectral resolution 
of A/ A A ~2000-4000. 

The GMOS data reduction pipeline was used to extract 
and wavelength-calibrate the spectra of each IFU element. 
The variations in fiber-to-fiber response were removed us- 
ing twilight flat-fields, and the wavelength calibration was 
achieved using a CuAr arc lamp. We also applied additional 
steps to improve the flattening and wavelength calibration 
using a series of custom idl routines. The wavelength cov- 
erage of the final data is typically 3900-5500A, but varies 
depending on exact instrumental setup. Flux calibration was 
carried out using observations of standard stars through the 
IFU using identical setups as the target observations. From 
the observations of the standard star, we note that the am- 
plitude of the point-spread function (PSF) varies by <0.05" 
across the wavelength range 4000-7000A and so the results 
in the following sections are insensitive to variations in PSF 
with wavelength (0.05" is ~l/4 of the size of a lenslet). 

Since the datacubes have a large wavelength coverage 
we correct for the parallactic angle in each cube by modeling 
a 3- hour observation running from —1 to +1 hour angle with 
the corresponding minimum and maximum airmasses. After 
building the datacube, we model and correct this aberration 
using a linear interpolation at each slice of the datacube 
along the wavelength axis. We note that the typical airmass 
of these observations was ~1.2 and the typical amplitude of 
the parallactic angle is ^0.3" across entire wavelength range 
of each cube. 



3 ANALYSIS 

3.1 Stellar Masses and Morphologies 

Before discussing the spatially resolved properties, we first 
make use of the existing imaging to estimate the stellar 
masses and morphologies of the galaxi es in our sam p le. To 
estimate the stellar masses, we follow iMcGee et aD <|201lh 
and use SED mod eling employing the Sloan ugriz- and K- 
band photometry dBalogh et al.lf2 005V Brie fl y, we generate 
galaxy templates using Bruzual fc Chariot] (l2003h models 
with a Salpeter initial mass function (jSalpeterl 1955h assum- 
ing a lower and upper IMF mass cutoff of 0.1 and 100 Mq 
re spectively. We follo w the galaxy parameter ranges used 
bv lSalim etaD l|2007t l. The model spectra include a range of 
galaxy age (0.1-13 Gyr), metallicity (0.005-2.5 Z ), star for- 
mation history and dust obscuration (Au=0-6 mags). In par- 
ticular, the exponentially declining star formation rates with 
superimposed bursts are randomly chosen with a uniform 
distribution, and we allow bursts that last some time ran- 
domly distributed in duration between 30 and 300 Myr. The 
strength of the bursts are also randomly chosen so that dur- 
ing the lifetime of the burst they produce between 0.03 and 

4 times the stellar mass the galaxy had at the onset of the 
burst. From these model star-formation histories, we gener- 
ate magnitudes by convolving the resulting spectra with the 
ugriz and A"-band filter response curves. To fit to the obser- 
vations, we search the entire parameter space and minimise 
the x 2 , but generate a probability distribution functio n for 
each parameter for each galaxy fsee lMcGee et af1l201ll for a 
detailed discussion). The resulting stellar masses (and their 
la errors) are given in Table 2. The median stellar mass of 
the galaxies in our sample is 8±2xl0 10 M Q . 

Since the K-band luminosity is a good tracer of the 
stellar mass (and is less sensitive to the recent star-formation 
history than the optical bands), we can also perform a simple 
check of the stellar mass using a canonical mass to light 
ratio. Bal ogh et al.l |2005l ') (see also Fig. [2]) show that the 
(u — g) and (r — i) colours of E+A galaxies are consistent 
with a truncated star-formation model in which a burst of 
5-15% is superimposed onto a old population, and use this 
to derive a canonical M/La"=0. 8+0.1. Applying this to our 
Jf-band magnitudes we derive M*=9±2xl0 10 Mq, which is 
consistent with the estimates from the more sophisticated 
modeling above. This stellar mass is also com parable to the 
averag e stellar mass of the parent sample from lBalogh et all 



(120051 ); M, =8±2xl0 10 M o . 



The SDSS imaging is also useful for measuring Sersic 
indices, asymmetries and radial colour gradients within the 
galaxies which, along with the distribution of H<5 equiva- 
lent widths, may provide a diagnostic of the physical mech- 
anism^) of E+A formation. For example, positive colour 
gradients (i.e. bluer in center) can arise when young stars 
are more concentrated than the old stellar population. In 
this case, the distribution of US equivalent widths are also 
likely to be compact. In contrast, a negative colour gradient, 
together with a positive H<5 equivalent width radial gradi- 
ent, may represent a galaxy in which the molecular clouds 
are not confined to the nuclear regions. 

From our sample of eleven galaxies, three show positive 
colour gradients (A(g — i)/Alogr >0.0), and the rest are 
eithe r flat or mildly negative (see also Fig. [71 and I Yang et al.l 
2008). The galaxies with bluer nuclei (positive color gradi- 
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Figure 3. Results from the IFU observations of the E+A galaxies in our sample. From top to bottom: 1. True colour BI K-h&i\d images 
of each of the galaxies constructed from Gemini/GMOS imaging (blue), SDSS i-band (green) and UKIRT/UFTI A'-band imaging (red). 
In each panel we overlay the GMOS IFU field of view. The arrows denote the directions of North and East. 2. Continuum images from the 
continuum above the 4000A break, derived by collapsing the datacube between 4350-4750A in the rest-frame, with dark grey indicating 
high intensity. The contours trace the continuum from the rest-frame 3650-3850A continuum emission (which should have a stronger 
contribution from young stars than the continuum emission above the 4000A break) . The solid bar in each panel denotes the seeing disk 
for the observations. 3. Gas phase emission line intensity (derived from either the [Oll]A3727 or [Oiii]A5007, as indicated). The contours 
start at 3<r and are incremented by la. 4. Two dimensional stellar velocity fields of the galaxies, measured from the Hi5 line. 5. Two 
dimensional stellar velocity dispersion maps. The dotted line marks the major kinematic axis in the velocity fields, where it can be clearly 
identified. 



ents) may evolve into the negative gradients typical in E/SOs 
if the central pa rts of these galaxies are metal-enhanced 
(|Yang et alj|2008h . Finally, we note that the median asym- 
metry and Sersic index of the sample is 0.03±0.01 and 
4.9±1.3 respectively, characteristic of bulge-dominated sys- 
tems. 

3.2 Spatial Distributions 

To investigate the spatial distribution of the gas and stars 
in the galaxies, first we extract narrow-band slices from the 
datacube around the emission and absorption lines of in- 
terest. For a given emission or absorption line, we fit and 
subtract the continuum using a region ±100 A from the cen- 
ter of the emission or absorption feature using a 3-<j clip in 
the fit to be sure that neighboring emission and absorption 
lines are omitted. We also extract continuum images from 
the datacube. For each image, we median filter each spectral 
pixel in the datacube between 4350-4750A (rest-frame) as 
a proxy for the "old" stellar population (above the 4000A 



break) and between 3650-3850A (rest-frame) as a proxy for 
the "young" stellar population and show these in the second 
row of Fig. [3] We caution that the continuum emission from 
above and below the 4000A break has contributions from 
both "old" and "young" stars. For example, for a 10% by 
mass starburst which has been truncated at 1 Gyr, for solar 
metallicity the contribution to the continuum emission from 
the burst to the total 3650-3850A and 4350-4750A contin- 
uum emission is ~60% and 40% respectively, but can vary 
from 40%- 7 0% in each band depending on starburst age and 
metallicity (|Marastonlll998l 'l. 

A better test of how the A-stars stars are related to 
the underlying stellar emission, (as well as the current star- 
formation activity) is to turn to the US absorption line. We 
estima te the equivalent w idth of H8 a t each spatial pix el fol- 
lowing [Baiogh^FaLl (|l999l ) (see also lGoto et al.l 120031 ). We 
first average each spatial pixel over a 0.6x0.6" region to ap- 
proximately match the seeing, and estimate the continuum 
flux using 3<r clipped linear interpolation between two wave- 
length windows placed at either side of the H<5 line (4030A 



6 Swinbank et al. 




-[0 II] intensity 



stellar velocity 



9 



Hps., 



E+A Jl 627+4800 




• 




* 





8 6 4 2 -2 -4 -6 - 
Continuum 




stellar velocity 



/s 20*0 




[0 II] 


in tens 














stellar velocity 



E+A J2307+1525 






• 










8 6 4 


2 0-2 


-4 -6 -8 



Continuum 

9 


jLO ill) intensity 

i . 








stellar velocity 

* 


* 

-40 km/s 40-. 


stellar a 

i 


20 km/s 150-.. 



-3-2-10 1 2 3 
Arc Seconds 



-1 1 
Arc Seconds 



to 4082A and 4122A to 4170A: lGoto et al.ll2003h . This 3x3 
averaging produces an effective seeing disk of 0.85", and 
this is taken into account in all of the following calcula- 
tions. The rest-frame equivalent width of the H<5 absorption 
line is then calculated by summing the ratio of the flux in 
each pixel of the spectrum, over the estimated continuum 
flux in that pixel based on our linear interpolation; its dis- 
tribution over the area of each galaxy is shown in Fig. [4] 
In this figure we also indicate the centroids of the nebu- 
lar emission and continuum from above the 4000A break, 
with a la error ellipse. To calculate these error ellipses we 
compute the number of photons per pixel with their as- 
sociated \fn uncertainties, and re-compute centroid with 
10 4 monte-carlo realisations. It is interesting to note that 
the regions of strongest H<5 absorption do not always coin- 
cide with regions with the strongest emission lines. Indeed, 
Fig. [4] shows that, five objects have H<5, nebular emission 
and continuum emission which are spatially co-located on 
~kpc scales (J0948,J1642,J1627,J1715,J2307), whilst three 
systems have significant offsets between nebular emission 
and US (J0906,J1242,J1348). Two galaxies have nebular 
emission and H<5 in agreement, but with very extended H<5 
which is offset from the strongest continuum (J0835,J1013). 
Thus, approximately half the sample are (relatively) sim- 
ple in their A-star/gas/continuum emission, whilst half 
are clearly more complex. On average, the centroid of 
the A-stars and the continuum emission are separated by 
ArHa- s tars=0.8±0.4 kpc, whilst the average offset between 



the Hr5 and nebular emission is Arns-gas— 1-0±0.2 kpc, but 
can be as large as as 1.7 kpc. This suggests that at least in 
some cases, the new OB stars are either displaced from the 
older A-stars, or that the extinction is not uniform. In either 
case this helps to reconcile the fact that both emission and 
strong absorption can be seen in the same, integrated galaxy 
spectrum. The two populations do not generally arise from 
exactly the same place. 

Next, we examine the spatial extent of the A-stars 
through the H<5 absorption. In the following calculations, 
we assume that both the continuum emission is smoothly 
distributed and subtract the amplitude of the point-spread 
function at the redshift of the galaxy in quadrature from 
the area's covered by the continuum and A-stars. As shown 
by the distributions in Fig. [4] the A-stars are typically dis- 
tributed over a large radius. We quantify this by computing 
the radius at which the H<5 equivalent width drops to half the 
peak value (this is extracted along the major kinematic axes 
within the galaxy which is discussed in § [3T3J . This quantity, 
rh(W H5) is compared with the observed g — i colour gra- 
dient within the central 2 kpc in the top-left panel of Fig. [7] 
Positive radial colour gradients indicate galaxies that are 
bluer in the center, and are measured in only three galax- 
ies (~30%). This is lower than the fraction of E+A galaxies 
wi th strong positive c o lour g radients (i.e. blue cores) found 
bv lYamauchi fc Goto! i|2005f ) (~67%), but still higher than 
the fr action of early type galax ies with positive colour gradi- 
ents l|Balcells fc Peletiedll994h . Within our sample, there is 
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Figure 4. The distribution of rest-frame H<5 equivalent widths across each of the galaxies in our sample. The contours are separated 
by A Wo (H5)=1A. The 'X' denotes the center of the galaxy as defined from the center of the continuum above the 4000A break, while 
the ellipses denote the center and la uncertainty in the nebular emission centroid (red) and continuum emission centroid (black). This 
figure highlights that six E+As in our sample have H<5 absorption, nebular emission and continuum emission which are spatially co- 
located on ~kpc scales (J0948, J1642, J1627, J1715, J2307), whilst three systems have significant offsets between nebular emission and H<5 
(J0906, J1242, J1348). Two galaxies have nebular emission and H<5 in agreement, but with very extended H<5 and nebular emission which 
is offset from the strongest continuum (J0835, J1013), suggesting that the A and OB stars formed in a different location than the older 
stars. At least in some cases, the new OB stars are either displaced from the older A-stars. This reconciles the fact that both emission 
and strong absorption can be seen in the same, integrated galaxy spectrum - the two populations do not generally arise from exactly 
the same place. 



no evidence for the correlation between colour gradient and 
spatial extent of the A-stars (as might have be en expected 
from merger or tidal models for E+A formation; iBekki et al.l 
l2005h . 

Another way to characterise the extent of the HW ab- 
sorption feature is to measure the area over which the equiv- 
alent width exceeds some threshold. We choose a thresh- 
old of 6A, and calculate the physical area A(H<5>6A) 
for each galaxy. The results are shown in Fig. [7] as a 
function of various kinematic parameters, discussed below. 
Note that in all but one galaxy (J0948+0230) the area 
with W (H<5)> 6A is significantly larger than the see- 
ing. The median A(H<5 > 6A)=4.0±2.1 kpc 2 , with a range of 
A(H(5 > 6A)=0.8-15.2 kpc 2 . In comparison to the continuum 
emission from above the 4000A break, the A-stars cover at 
least ~10% of the area, with a median of 33±5%. Such wide- 
spread A -star popula t ions h ave also been found in previous 
samples: iGoto et al.l l|200ct) show that the A-stars in two 
SDSS selected E+A galaxies are di stributed over ^10 kpc 2 , 
while the seven E+A galaxies from iPracv et all ((2009) have 
A-stars that are typically distributed over 5-15 kpc 2 . 

In summary, we note two important observa- 
tions. First, the Hr5 is widespread, typically over 
A(H(5 > 6A)=4. 0+2.1 kpc 2 and thus the unusual spectrum 
characteristic of E+A galaxies is a property of the galaxy as 
a whole. It is not due to a heterogeneous mixture of popu- 
lations, such as a nuclear starburst with a very high equiv- 
alent width (>10A) on top of an otherwise normal galaxy. 
Second, in approximately half of our sample, the A-stars, 
nebular emission and continuum are not co-located, sug- 
gesting that the newest stars are often forming in a different 
place than those that formed glGyr ago, and that recent 



star-formation has occurred in regions distinct from the con- 
tinuum emission (the older stellar populations) . 

3.3 Kinematics 

To obtain the stellar velocity structure within the galaxies 
we first construct a continuum image of the galaxy from 
the datacube collapsed between 4350-4750A, and bin the 
data to a constant sign al-to-noise ratio per bin u sing the 
Voronoi binning method (|Cappellari fc Copidl2003r i . We de- 
mand a signal-to-noise of 30 in continuum, and then cross 
correlate each of the indi vidual spatial regions with a li- 
brary of template spectra l|Vazdekis fc Arimotd ll999l. The 
cross correlation is performed using the penalized pixel fit- 
ting method (ppxf; ICappellari fc Emsellemll2004l ). To max- 
imise signal-to-noise, to derive the velocity field of the young 
stellar population we use the wavelength range ~3800-4290, 
which covers Ca H&K and H<5. 

As Fig. [3] shows, all eleven galaxies have velocity gra- 
dients in their dynamics, with amplitudes ranging from 30- 
200 km s -1 . At least eight have dynamics that appear to be 
consistent with rotation; i.e., the peak of the line of sight 
velocity dispersion is at the dynamical center, and the ve- 
locity field is characteristically bisymmetric, showing a so- 
called "spider" diagram. For these galaxies, we model the 
velocity field in order to estimate the inclination and true 
rotational velocity. To achieve this, we use a simple ap- 
proach and model the velocity field with an arctan func- 
tion to describe the shape of the rotation curve such that 
v(r) = v c a,rctan(r/r t ), where v c is the asymptotic rotational 
velocity and r t is the effective ra dius at which the rotation 
curve turns over l|Courteaull 19971 ). We then construct a two- 
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Figure 5. The dynamical measurement of i>sin(i)/<7 as a func- 
tion of absolute i?-band magnit ude for E+A galaxi es from our 
sample and the recent work by iPracv et al.l d2009li The solid 
black points denote the average values in bins of Am=0.5 and 
the dashed line shows a linear fit, indicating a weak trend such 
that fainter galaxies tend to have higher v sin(i)/a. The plot also 
inclu des a comparison to t he early-type and elliptical galaxies 
from I Emsellem et al.l J2007T) (we have colour corrected their sam- 
ple from Ms to using a (B — R)=2.0 appropriate for an 
early type galaxy at z=0.05). This sample have a lower (v/a), 
with vsin(i)/a (0.30±0.03), but show a comparably weak trend 
of with increasing absolute magnitude. Galaxies which are weak 
radio sources (and thus may harbor AGN are shown as open sym- 
bols); do not stand out from the rest of the sample. 



dimensional kinematic model with six free parameters (v c , 
r t , [x/y] center, position angle (PA) and inclination) and 
use a genetic algorithm with 10 random initial values and 
a scale factor of 0.95 to search for a best fit. We demand >30 
generations are performed before testing for convergence to 
a solution, throwing out (a maximum of) 10% of the solu- 
tions in each generation. The best fit kinematic maps are 
shown as contours in Fig. [3] and the best-fit inclination and 
rotation speed are given in Table 2. Using the velocity maps 
and model fits, we identify the major kinematic axis wher- 
ever possible and extract the one dimensional rotation curve 
and velocity dispersion profile along and show these in Fig. [5] 
From each of the two-dimensional maps, we also extract the 
distribution of H<5 equivalent widths, both in radial bins and 
along the kinematic cross section. 

A measure of whether the dynamics are dominated by 
ordered- or random- motion can be obtained by measur- 
ing the ratio of rotational- velocity to line of sight velocity- 
dispersion, v sin(i)/cr. We derive a median luminosity- 
weighted v sin(i)/<r=0.8±0.4 for the whole sample with a 
range of 0.2-5.5. However, Fig. [5] shows that there is a 
weak correlation between absolute magnitude and v/a, such 
that the lower luminosity systems have an increasing ro- 
tational support, although only four of the eleven galaxies 
have v/a > 1, uncorrected for inclination. Recen t studies 
of early type galaxies from lEmsellem et al.l (|2007l ) identify 
a comparably weak trend of v/a with magnitude (although 
we note that their sample has a slightly lower average v/a 
wi th an average v/cr=0. 30±0,03). Indeed, the early types 
in lEmsellem et ail (|2007r) have Mk=-2 1.8+0.3 compared to 



M_r=-21.5±0.5 for the E+As, and with stellar masses of 
1.0±0.3x 10 11 Mq, which is also consistent with our sample. 
Thus, these early- types and ellipticals from lEmsellem et ail 
(2007) also have many of the properties expected for the 
des cendant popul a tion o f E+A galaxies. 

iBender et ail (| 19941 ) also find a correlation between 
magnitude and v/a for lumin ous elliptical galaxies. Al- 
though the IBender et ail (1 19941 ) sample are dominated by 
brighter galaxies than the E+A's studied here (typically 
~ 1.3 mags brighter), those galaxies with high v/a (~1 and 
so comparable to the E+A's) are increasing dominated by 
lower mass, lower luminosity Sb-SOs. 

A better test of the dynamical state of galaxies with 
low v/a can be made by measuring the luminosity- weighted 
ste llar angular momentum per unit mass, Xr as developed 
by lEmsellem et ail (120071 ). This combines the velocity am- 
plitude and line of sight velocity dispersion via: 
Xr —< R\v\ > / < R\/v 2 + a 2 >, where R, v and a are 
the radius, rotational velocity and velocity dispersion from 
the kinematic center. Xr is designed to quantitatively distin- 
guish galaxies that have similar v /a values but very different 
velocity structures by identifying those g alaxies with sig- 
nifican t angular momentum per unit mass. lEmsellem et ail 
(|2007l ) define "fast" and "slow" rotators as Xr >0.10 and 
Xr < 0.10 respectively. Classically, "fast" rotators are domi- 
nated by low mass spheroids (with a large v/a) whilst "slow" 
rotators are dominated by ellipticals, albeit with a large 
range in M_r. We report Xr for each galaxy in Table 3. Our 
sample has a median A jj=0. 35+0.1, and all of the galaxies in 
our sample have Xr > 0.1. When including the E+A's from 
IPracv et ail l|2009l) . 18/19 galaxies with well resolved dy- 
namical maps have Xr > 0.1. Taken with the measurements 
oiv/a, Xr and stellar masses, the E+As in our sample have 
properties comparable to those measured in local, represen- 
tative elliptical galaxies and SOs where similar measurements 
have been made. 

To investigate how the dynamics and A-star popula- 
tions are related, we correlate the kinematics as measured 
from v/a and Xr as a function of the A-star spatial ex- 
tent and colour gradients. Fig. [7] shows that the galaxies 
which have the highest v/a also have the most compact 
A-star distributions. Indeed, the average v/a for the galax- 
ies where A(H<5 > 6A)>4kpc 2 is 0.5+0.2, compared with 
v/a—1. 1+0.5 for the remainder. We note that we would ex- 
pect the opposite correlation between the spatial extent of 
the A-stars and v/a if the measurements were being driven 
by signal-to-noise effects alone (low v/a values would be 
expected for the most compact sources where it is more dif- 
ficult to identify large velocity gradients). However, mea- 
suring the spatial extent of H<5 absorption is sensitive to 
detecting the continuum, so we also normalise the area with 
A(H<5 >6A) by the spatial extent of the continuum emission, 
A(gal), as measured from the continuum. For the continuum 
measurement we use the flux redward of the 4000A break as 
measured from our IFU observations, although we note that 
the results are quantitatively similar when using the r-band 
half light radius, or R e , measured from SDSS imaging. In all 
panels, we find that galaxies with high v/a show the most 
compact A-star populations. 

Of course, the spatial distribution of an equivalent 
width depends on its intrinsic spatial distribution as well 
as the overall distribution of the continuum light (or galaxy 
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Figure 6. Extracted, one dimensional equivalent width, velocity and line of sight velocity dispersion profiles are shown for each of the 
E+A galaxies in our sample. The profiles are extracted along the major axis as indicated in Fig. [3] Left: One dimensional extent of 
equivalent widths. Middle: One dimensional rotation curve of the galaxy with the best fit rotational model (where appropriate) overlaid 
as contours Right: One dimensional distribution of line widths. The solid bar in the left hand panel denotes the amplitude of the seeing 
disk, converted to physical scale at the redshift of each galaxy. 



surfac e brightness profi le) after convolving with the seeing 
(e.g. IPracv et al.l 120051 ). In particular, for a galaxy with a 
centrally concentrated equivalent width profile, the steeper 
the continuum surface brightness profile, the more the stellar 
population in the galaxy center will contamina te the outer 
regio ns after convolution with the seeing (see IPracv et al.l 
|2010| for a detailed discussion). In principle this could give 
a correlation between v/a and the extent of the A-star 
light distribution. Whilst full modelling to recover the in- 
trinsic equivalent width distribution is beyond the scope 
of this paper, we can perform a simple test of whether 
this effect is likely to significantly contribute our results. 
We therefore construct 10 4 mock GMOS datacubes with 
variable amplitude velocity fields, continuum emission- and 
H<5 equivalent width- profiles (which broadly bracket the 



range of our observations) and convolve these with the see- 
ing. We note that in these simulations, the extent of the 
continuum emission and H<5 profiles are independent. We 
measure the "true" v/a and ratio of A(H8 >6)/A(gal) as 
well as the same quantities from the recovered cubes after 
convolution with the seeing and surface brightness limits 
appropriate for our observations. Allowing the spatial ex- 
tent of the A-stars to vary from "unresolved" to 3x the 
continuum emission, we find that the change in slope of 
Alog(v/a)/A(A(HS >6)/A(galaxy)) is -0.04 to -0.2 depend- 
ing on choice of parameters. The gradient of the correla- 
tion seen in Fig.[7]is Alog(v/a)/A(A(HS >6)/A(galaxy))=- 
2.2T0.3. Thus the underlying A-star distribution, contin- 
uum surface brightness profile and seeing may contaminate 
the correlation between the extent of the A-stars and dy- 
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Figure 7. The spatial extent of the A-stars within the E+A galaxies in our sample as a function of colour gradient and kinematics. 
Top Left: The extent of the A-stars (defined as the radius at which the H<5 equivalent width drops to half the peak value along the 
major kinematic axis) as a function of the g — i colour gradient within the central 2 kpc of the galaxy. Positive radial colour gradients 
denote light that is bluer in the center. The dashed line denotes a colour gradient of zero, and we use this to differentiate galaxies with 
positive- (black filled circles) from negative- colour gradients (red filled squares). In all panels we identify the E+As which have weak 
radio emission (L1.4 ^10 22 W/Hz; highlighted with an open circle). Top Right: The spatial extent of the A-star s as a func t ion o f the 
dynamics as measured from (v sin(i)/a) for the E+A galaxies in our sample. We also include the six E+As from lPracv et al,l 1120091 ). To 
provide a quantitative measure of the spatial extent of the A-stars, we define A(H5 > 6A) as the area over which the H<5 is stronger than 
6A. The galaxies with higher v sin(i)/o tentatively show signs of having more compact A-star distributions compared to galaxies which 
are pressure supported; the average tisin(j)/<r for the galaxies where A(H(5>6A)>4 kpc 2 is 0.5+0.2, compared with v sm(j)/<r=l.l±0.5 
for the others. Bottom Left: Since the measurement of the spatial extent of H<S is sensitive to detecting the continuum, we normalise the 
area, A(H5), by the spatial extent of the galaxy, as measured from the continuum above the 4000A break. We note that the results are 
insensitive to whether we use the SDSS r-band image to measure the half light radius, or R e . This plot shows that the galaxies with the 
most compact H5 distribution tend to be those which have dynamics consistent with rotation, while the galaxies with A-stars which are 
most wide-spread are those with disp ersion dominated dyna mics. Bottom Right: The spatial extent of the A-star population normalised 
to the continuum as a function of A^ iEmsellem et al1l2007tl . High and low \r correspond to 'fast' and 'slow' rotators respectively. This 
plot shows that the fastest rotators tend to have the most compact A-star spatial distribution. We note that although our sample is 
small, the galaxies with weak AGN do not stand out in terms of their A-star extent or dynamics from the rest of the sample in any of 
the panels. 



namics seen in Fig. but it seems unlikely that this effect 
alone can account for the correlation we see. We also note 
that in general, rotation-dominated galaxies have shallower 
surfa ce brightness pr ofiles than pressure-supported galaxies 
(e.g. ICourteau|[T997h . This means there is a physical corre- 
lation between the two which may lead to an indirect cor- 
relation between v/a and A(H<5)/A(gal). However, clearly, 
observations of a well matched control sample of non E+A's 
around the H<5 are required to examine whether the corre- 
lation in Fig. [7] is a result of the processes which drive the 
E+A's. Nevertheless, for our sample the 'slow-rotators' tend 



to have the most wide-spread A-star populations, while the 
'fast-rotators' have the more compact A-star distributions. 



4 DISCUSSION 

There are at least three important unresolved problems re- 
garding E+A galaxies: 1) from where in the galaxy does 
the unusual spectrum originate? 2) what star formation his- 
tory leads to the spectral characteristics of these galaxies; 3) 
what triggers the recent change in star formation history? To 



Integral Field Spectroscopy of E+A Galaxies 11 



Table 2: Properties of the E+A Galaxies in Our Sample 





Equivalent 
[On] 
(A) 


Width 
H<5 

(A) 


B/T 


Sersic 
index 
(") 


(X10 10 
Mg) 


[On] Flux 
( x 10 _ 16 
crg/ S /cm 2 ) 


SFR([On]) 
(Mq /yr 




(deg) 


(km/s) 




r(H5) 
(kpc) 


(kpc'' ' 


J0835+4239 


2.9±0.2 


6.9±0.3 


0.48 


7±1 


«±?:Z 


5.0±0.5 


0.14±0.07 


0.4±0. 


1 


51±11 


55±15 


0.11±0.04 


2.4±0.2 


15.2 


J0906 + 5221* 


1.5±0.3 


6.3±0.4 


0.69 


4±1 


a . +0.5 
23.4±6.0 


2.3±1.0 


0.08±0.04 


1.0±0. 


1 


26±5 


75±12 


0.27±0.04 


2.0±0.2 


11.4 


J0938 + 0001 


3.9±0.4 


5.0±0.3 


0.47 


3.5±1.0 


5.1±0.5 


0.15±0.04 


5.5±0. 


1 


57±6 


202±12 


0.67±0.05 


1.2±0.5 


>0.8 


J0948 + 0230 


2.4±0.5 


6.4±0.2 


0.29 


2.2±0.3 


6.6+1- 


3.8±0.9 


0.05±0.03 


1.8±0. 


1 


54±7 


200±10 


0.53±0.03 


1.2±0.4 


>0.8 


J1013 + 0116 


5.2±0.3 


B.8±0.5 


0.46 


>7 


±11 


12.1±0.2 


0.48±0.12 


0.2±0. 


1 


42±10 


152±14 


0.27±0.04 


2.8±0.3 


14.4 


J1242 + 0237* 


2.0±0.4 


6.0±0.2 


o. to 


3.8±0.5 


i..+o : « 


2.7±0.6 


0.06±0.Q3 


0.5±0. 


1 


21±5 


92±6 


0.20±0.03 


1.0±0.2 


2.1 


J1348 + 0204 


2.1±0.5 


6.0±0.3 


0.67 


6.2±0.4 


...+*;■ 


7.1±0.8 


0.11±0.03 


0.4±0. 


1 


49±7 


63±8 


0.10±0.04 


1.6±0.2 


2.2 


J1627+4800 


0.8±0.4 


5.2±0.5 


0.73 


>7 




<2.1 


<0.13 


2.2±0. 


1 


37±5 


229±12 


0.52±0.05 


2.1±0.5 


4.4 


J1642+4153 


1.2±0.5 


6.5±0.5 


0.64 


3.5±0.6 


6.0+;- 


4.5±1.0 


0.08±0.03 


1.2±0. 


1 


51±6 


201±10 


0.41±0.04 


1.5±0.3 


3.2 


J1715 + 5822 


2.2±0.4 


6.5±0.8 


0.48 


4.5±1.0 


™t?:Z 


2.8±0.5 


0.17±0.05 


0.4±0. 


1 


49±5 


90±5 


0.46±0.03 


2.7±0.3 


6.5 


J2307+1525 


1.7±0.3 


6.3±0.4 


0.45 


>7 




3.0±0.5 


0.05±0.02 


0.7±0. 


1 






0.22±0.03 


1.4±0.2 


2.8 



Table 2. Notes: The two E+A galaxies labeled with a * (J0906+5221 and J1242+0237) are both in cluster environments, n denotes the 
Sersic index measured from the SDSS gri photometr y. The estimate o f the star- formation rate (SFR) is calculated using the [On]A3727 
emission line flux and assuming the calibration from lKennicut j il998h . 



address these issues, we have performed three dimensional 
spectroscopy of a sample of eleven massive E+A galaxies 
selected from the SDSS for their unusually strong H5 equiv- 
alent widths but weak [On] emission, suggesting that star 
formation in these galaxies was recently truncated, possibly 
following a starburst. 

Our sample was selected to span a range of morphol- 
ogy and environment and hence is diverse by design. More- 
over, with a limited sample size of evelen galaxies it is not 
possible to answer any of these questions definitively, espe- 
cially since there may be more than one way to form an 
E+A galaxy. Nevertheless, there are several important re- 
sults that are generic within our sample. First, we note that 
we see no strong correlation of the galaxy dynamics with the 
K— band morphology, or the environment of the galaxy with 
their kinematics or morphology of the nebular emission or 
A-stars: the two E+A galaxies in dense environments do not 
appear to stand out from the rest of the sample in any of 
their properties, so we discuss the properties of the ensemble 
below. 

To investigate the likely age of the starbursts and hence 
star-formation history, we estimate the age of the A-star 
populations by comparing the average gala xy spectrum fo r 
each galaxy with the spectral library of lJacobv et al.l i|l984h . 
We use the strength of the Ca K 3 933A absor ption line 
(which is especially sensitive to age; lRosel ll985). together 
with the equivalent width of the Ca H 3969A (and He) ab- 
sorption lines. For all of the E+A galaxies in our sample, 
the spectra best resemble that of an A5+1 star, with an ef- 
fective temperat ure of 8160 K. Inte rpolating the theoretical 
isochrones from iBertelli et al.l (|l994T ) , this suggests a lumi- 
nosity weighted age of 0.8+0.1 Gyr (for the MS turn off at 
A8) although we caution that stellar mixes and metallicity 
make this estimate uncertain. Nevertheless, this is also sup- 
ported by the broad-band colours which are well matched 
by a model in which an instantaneous starburst (which ac- 
counts for 10% by mass), is superposed upon the continuum 
above the 4000 A break (Fig. [5]). 

Thus, the colours and line strengths we observe can 
best be reproduced with ~ 10 per ce nt (by mass) star- 
bursts on top of a n old population (e.g. iBalogh et al.ll2005l ; 
IShiovaet aill2004r i. It is also evident that any residual star 



formation in the galaxy is negligible: the median stellar mass 
of the E+A galaxies in our sample is 8+2 x 10 10 Mq , suggest- 
ing an average burst mass of Mb urs t ~1.0x10 10 Mq. The 
median star-formation rate of galaxies in our sample (in- 
cluding limits) is SFR([On])=0.13±0.03M Q yr" 1 , less than 
2% of the past average star-formation rate. Thus, the pres- 
ence of even moderately strong H<5 absorption (> 3 A or so) 
indicates that the star-formation rate at truncation was at 
least comparable to the past average rate. Integrating the 
current star-formation rate over 1 Gyr we derive a mass of 
~lx 10 8 Mq, or only ~ 0.1 per cent of the total galaxy mass. 
Since the mass involved in the burst is at least ~10% of the 
mass of the stellar mass, thus suggests that the main pro- 
genitors of our E+A galaxies had significantly higher star- 
formation rates than seen today, and were most likely blue- 
cloud, star-forming systems. 

The A-stars are widely distributed, with an average 
A(H5 > 6A)=4.0±2.0 k pc 2 , and a range of A(H <S >6A)=0.8- 
15.2 kpc 2 (see a l so ISwinbank et al.l 120051 : iGoto et al.l 
120081 : iPracv et all I2009T ). On average, the strongest H<5 
(W (tL5)>6A) covers 33+5% of the area defined by the 
detected continuum above the 4000A break. It is also in- 
teresting that none of our galaxies show regions with very 
strong absorption, W o (HW)>10A, at least on ~kpc scales. 
Such strong absorption would be an unambiguous indica- 
tion of a significant enhancement of the star formation rate 
in the past g 1 Gyr; weaker lines can be formed via trun- 
cation of star form ation in normal, gas-rich galaxies (e.g. 
IBalogh et ai] |l999). Moreover, it is also clear that in ap- 
proximately half of the sample, the A-stars (which give rise 
to the absorption) and nebular emission are not co-located, 
suggesting that the newest stars are often forming in a differ- 
ent place than those that formed ;£ lGyr ago. Thus, together 
with the burst mass to stellar mass ratio, this suggests that 
there needs to have been a recent, dramatic change in star 
formation history, coordinated over kiloparsec scales in these 
galaxies. 

The galaxies display a range of kinematic properties: all 
eleven have velocity gradients in their stellar kinematics and 
we derive a median luminosity-weighted v sin(z)/er=0.8±0.4 
for the whole sample (with a range of 0.2-5.5). Fig.[5]shows 
that there is a weak correlation between absolute magnitude 
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and v/a, such that the lower luminosity systems have an 
increasing rotational support. The median Xr for the E+A's 
in our sample is Aj?=0. 35+0.1, and all of the E+A galaxies 
studied here have Xr >0.1 which is consistent with that 
found for a comparably bright sample of elliptical galaxies 
and passive SO' s which are also well m atched in their average 
stellar masses jEmsellem et alj|2007ft . 

Although we have shown that the E+A phase is a 
galaxy-wide phenomenon, conclusively addressing what trig- 
gers the E+A phase is more difficult, especially since there 
are a number of processes which could be responsible at once 
(e.g. strong interactions or mergers, possibly proceeded by 
AGN activity). However, we can search for global trends in 
the A-star population with other galaxy properties in order 
to search for the major contributing factor. We find that dy- 
namics are correlated with the spatial extent of the A-stars 
such that the 'slow-rotators' have the most wide-spread A- 
star populations, while the 'fast-rotators' have compact A- 
star distributions. Since the E+A galaxies in our sample do 
not preferentially lie in dense/cluster environments, it is un- 
likely that any truncatation of star-formation was due to 
ram pressure stripping of the gas. 

It is interesting to note that theoretical models which 
aim to trace the formation and evolution of early type galax- 
ies have shown that there are two primary factors which 
determine the the structure of early-type galaxies. More 
anisotropic and slowly rotating galaxies result from pre- 
dominantly collisionless major mergers, while faster rotating 
galaxies are produced by more gas -rich mergers, where dis- 
sipat i on plays an important role |Kprmendv fc Diorgovskil 
1 19891 ; iBender et aflll992l ; iFaber et alj 1997ft . The second ma- 
jor factor which drives the galaxy structure is the mass frac- 
tion of the merger components. Unequal-mass galaxy merg- 
ers tend to produce more rapidly rotat ing galaxies than 
mergers which compa rable mass ratios (|Naab et al.l 1 19991 ; 
iNaab fc Burkerj|2003ft . Whichever process was responsible, 
the encounter has left residual rotational motion in at least 
90% of our sample, with A-stars widely distributed within 
the ISM. Qualitatively, our results are consistent with mod- 
els of unequal mass, gas rich mergers. 

Finally, it is interesting to examine whether activation 
of a cent ral black hole cou ld be related to the end of star for- 
mation (|Yan et alj |2006a). As described in § 12.21 four of our 
galaxies are detected in FIRST with 1.4 GHz fluxes in the 
range 0.5-2.2 mjy, corresponding to luminosities of Li.4=8- 
15xl0 21 WHz~\ All four galaxies with FIRST detections 
show [Oin] and (weaker) [On] emission. For these galaxies, 
the centroid of the nebular emission appears offset from the 
continuum by 0.4-1 kpc, but this is consistent with the rest 
of the sample. Thus, if there is an AGN in these galaxies, 
it does not appear to be the sole contributor to the nebular 
emission. Moreover, if the radio emission were due to star 
formation, then t he instantaneous s tar-formation rate would 
be 6-1 1M yr" 1 (|Helou et al.ll 19851 ) a factor ~60 times that 
inferred from the nebular emission lines. This could be pos- 
sible if the galaxies are exceptionally dusty, although this 
seems unlikely given the blue (r — i) colours; we would ex- 
pect dusty, star-forming galaxies to lie on the up per edge of 
the blue cloud in Figure [H with (r - i) > 0.4 (jWolf et all 
120051 ; iBalogh et a"Hl2009ft . Thus it is most likely that the 
radio emission arises from an AGN. In all the key figures 
we have highlighted the E+A galaxies in our sample that 



have Li.4 ^10 22 W/Hz. In all cases their properties appear 
to be consistent with those of the rest of the sample. Thus, 
although we only have small number statistics, our results 
suggest that AGN activity in E+A galaxies does not play a 
dramatic role in defining the spatial or kinematic properties 
of the A-stars, and/or that its effects are short-lived. 



5 CONCLUSIONS 

We have investigated the spectro-photometric properties 
of a sample of eleven E+A galaxies in the redshift range 
0.12 < z < 0.07 in order to provide constraints on the phys- 
ical processes which cause the unusual post-starburst signa- 
tures. Although our sample selection was hetrogeneous (by 
design we selected E+A's with a range of morphologies and 
environments), there are several important results which can 
be summarized as: 

• The strongest H<5 (W >6A) tends to widely dis- 
tributed within the galaxies, on average covering ~ 33% 
of the galaxy image and extending over areas of 1-15 kpc 2 . 
This suggests that the characteristic E+A signature is a 
property of the galaxy as a whole and not due to a hetero- 
geneous mixture of stellar populations. 

• In approximately half of the sample, the A-stars, 
nebular emission and continuum are not co-located. The 
average offset between the H5 and nebular emission is 
ArH<5- g as = 1.0±0.2 kpc, but can be as large as as 1.7 kpc. 
This suggests that at least in some cases, the new OB stars 
are either displaced from the older A-stars, or that the ex- 
tinction is not uniform. In either case this helps to reconcile 
the fact that both emission and strong absorption can be 
seen in the same, integrated galaxy spectrum. The two pop- 
ulations do not generally arise from exactly the same place. 

• The colours and line strengths of these E+A galaxies 
are consistent with a ~ 10 per cent (by mass) starburst on 
top of an old population suggesting an average burst mass 
of Mburst ~1O 1O M0. The current (residual) star-formation 
is unable to account for the burst mass. This indicates that 
the main progenitors of most of our E+A galaxies were blue- 
cloud, star-forming galaxies in which there was a recent, 
dramatic change in star formation history, coordinated over 
kiloparsec scales. 

• The kinematics show that some level of rotation in 
the A-star population is common, and the dynamics (v/a 
and Xr) are consistent with those measu red in compara- 
bly br ight early types recently studied by Em sellem et al.l 
(2007), which are also well matched in their average stellar 
masses. Whichever process was responsible for causing the 
starburst and subsequent truncation, the encounter has left 
residual rotational motion in at least 90% of the sample, 
with A-stars widely distributed within the ISM. 

• We find that the A-star dynamics are correlated 
with their spatial extent such that the 'slow-rotators' have 
the most wide-spread A-star populations, while the 'fast- 
rotators' have compact A-star distributions. 

• We also find that the fraction of galaxies which have 
radio emission suggestive of low luminosity AGN is 20-40%, 
a factor ~8x higher than expected given their stellar masses, 
indicating a high AGN fraction of E+A galaxies. However, 
although our sample is limited to only eleven galaxies, the 
kinematics and spatial distribution of the stars in the ra- 
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dio detected sub-sample do not stand out from the radio- 
undetected systems, suggesting that AGN feedback in E+A 
galaxies does not play a dramatic role in defining their prop- 
erties, and/or that its effects are short. 

Whichever process was responsible for the E+A sig- 
natures rotational motion is seen in the majority of our 
sample, and their dynamics, stellar masses and luminosities 
are consistent with those expected for the progenitor pop- 
ulation of (at least a subset of) representative, early type 
galaxies. Overall, these observations provide new and de- 
tailed constraints on the kinematics and spatial distribution 
of the A-stars and gas in E+A galaxies, and hence con- 
straints on galaxy formation models which aim to test the 
variety of physical mechanisms which trace the evolutionary 
sequence linking gas-rich, star-forming galaxies to quiescent 
spheroids. 
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